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The masses and dispersions of light hadrons are calculated in lattice QCD using 
an 0(a 2 ) tadpole-improved gluon action and the D^34 action, an 0(a 2 ) tadpole- 
improved next-nearest-neighbor fermion action originally proposed by Hamber and 
Wu. Two lattices of constant volume with lattice spacings of approximately 0.40 fm 
and 0.24 fm are considered. The results reveal some scaling violations at the coarser 
lattice spacing on the order of 5%. At the finer lattice spacing, the calculated N/p 
mass ratio reproduces state-of-the-art results using unimproved actions. Good 
dispersion and rotational invariance up to momenta of pa ~ 1 are also found. The 
relative merit of alternative choices for improvement operators is assessed through 
close comparisons with other plaquette-based tadpole-improved actions. 



1 Introduction 

Lattice discretization of the continuum QCD action introduces errors at finite 
lattice spacing a. The standard Wilson gauge action has 0(a 2 ) errors and 
the standard Wilson fermion action has 0(a) errors. Simulations using these 
actions have shown that lattice spacings of 0.1 fm or less and lattice volumes 
of 24 4 or larger are needed in order to hold^ystematic errors to the 10% level. 
Such simulations are major undertakingsErEI and require enormous computing 
power to extract even the most basic of hadronic observables, the hadron 
masses. 

During the past few^ears, considerable efforts have been devoted to im- 
proving lattice actions era. The idea is to reduce or remove the discretization 
errors from the actions so that they have better continuum-like behavior. At 
the same time, errors due to the lattice regularization are accounted for through 
the renormalization of the coefficients multiplying the improvement operators. 
The hope is the use of improved actions will allow one to simulate efficiently 
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and accurately on coarse lattices, such that computer resources may be redi- 
rected to the simulation of QCD rather than quenched QCD. Moreover, one 
may turn the focus of investigation towards quantities of experimental interest. 

In the pure gauge sector, the 0(a 2 ) tadpole-improved action □ leads to 
dramatic improvement in the static potential and glueball massestl up to lattice 
spacings of 0.4 fm. In the light quark sector, hadron spectroscopy has been 
investigated, with a variety of improved actions_.including the 0(ftLimproved 
SW action 0, the 0(a 2 )-improved_p2M_action El and its variants Em, and the 
D%34 action of Hamber and WuE3'UE3E] considered here. 

The D%34 action is an 0(a 2 ) next-nearest- neighbor fermion action with 
tadpole-improved estimates of the coupling renormalizations. This action is 
selected primarily due to its simplicity. The cost of simulating it is about a 
factor of two as compared to standard Wilson fermions. Our goal is to study 
its feasibility as an alternative action to SW which has the clover term, or to 
D234 which has both next-nearest-neighbor couplings and the clover term. In 
particular, we examine dispersion relations and test the rotational symmetry 
of both the gauge and the fermion actions. Hadron mass ratios are calculated 
for a wide variety of hadrons including hyperons. To explore scaling violations, 
we consider two coarse lattices of approximately fixed physical volume: 6 3 x 12 
at a lattice spacing of 0.40 fm and 10 3 x 16 at 0.24 fm. 

2 Improved Lattice Actions 

The improved gauge action employed in this investigation is given byi: 

So = /?E l Rc Tr (! - U pl) -iEi Re ^ - U rt) ■ 
pi rt 

The second term removes the 0(a 2 ) errors at tree level. Perturbative correc- 
tions are estimated^ to be of the order of 2-3%. Here, U T i denotes the rect- 
angular 1x2 plaquettes. uq is the tadpole factor that largely corrects for the 

large quantum renormalization of the links U^(x) = exp(z g J^ +a>1 A(y) ■ dy). 
In this calculation we use the mean plaquette u = (Re Tr(J7 pl )) 1 / 4 /3 to esti- 
mate uo, and will focus our evaluation of lattice action improvement on other 
plaquette-based improved actions, uq is determined self-consistently in the 
simulation. 

The improvement program of Sheikoleslami and Wohlert (SW)Ej provides 
a systematic approach to the improvement of lattice fermion actions. However, 
the on-shell improvement program leaves some freedom in the relative values of 
the coefficients of the improvement operators. In this investigation, we consider 
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Figure 1: Free dispersion relations for zero mass Dx34, Wilson, and continuum fermions. 
Momentum p is along the (1,1,0) direction. Beyond the D%34 branch point, the real part of 
the two conjugate roots is shown. 

a specific case of the general class of D234 actions i in which the improvement 
parameters are tuned to remove the second-order chiral-symmetry-breaking 
Wilson term at tree level. This fermion action may be written 

M Dx34 = m q + 7 • V + - (-a 2 V M A M + b a 3 A* ) (2) 

where 

V M t/,(x) = -!— [U^x) ^{x + fJ,) — Ul(x - (j,) i,(x - M )] , (3) 
2a uo F 

and 

A M ip(x) = [U^x) ij)(x + n) + UUx - /i) ijj{x -n)-2u i/)(x)] . (4) 

The second-order term of the D234 action, ^ M A M + a ■ F/2, breaks chiral 
symmetry and does not appear in the D%34 action. The fourth-order term 
of (^|) breaks chiral symmetry and provides for the removal of the fermion 
doublers. The D%34 action is free of both 0(a) and 0(a 2 ) errors at tree 
level. Explicit evaluation of (g) combined with a Wilson- fermipnpatade field 
renormalization, discloses the fermion action of Hamber and WuE3EJIl3. 

Free dispersion relations can be obtained by locating the poles in the 
fermion propagator. Fig. [j] shows free dispersion relations for massless con- 
tinuum, Wilson, and D^34 fermions. It is clear that D%34 fermions follow 
the continuum more closely than Wilson fermions. Note that there exists an 
unphysical high energy doubler (or ghost) in the D%34 action. It is very sim- 
ilar to that for the D234 action □ and is a general feature of fermions with 
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Figure 2: Static potential from Wilson loops. The empty circles are for j3 = 6.25, the solid 
circles for /3 = 7.0. The lines are best fits. The statistical errors are from 200 configurations 
in both cases. 



next-nearest-neighbor couplings. The doubler can be 'pushed away' from the 
low momentum region by various techniques, such as tuning the value of 6, 
or using an anisotropic lattice LJ. We simulate with 6=1, which gives good 
dispersion to pa ~ 1 as illustrated in Fig. |l|. 

3 Lattice Simulations 

3.1 Methods and Parameters 

Quenched gauge configurations are generated using the Cabibbo-Marinari pseudo- 
heat-bath methodEj. Periodic boundary conditions are used in all directions for 
the gauge field and in spatial directions for the fermion field. Dirichlct bound- 
ary conditions are used for the fermion field in the time direction. Configu- 
rations separated by 300 sweeps are selected after 4000 thermalization sweeps 
from a cold start. 

The static potential is calculated and the string tension a is extracted 
from the ansatz V(R) = Vq + oR — E/R where Vq and E are constants. 
Fig. H shows our results. Good rotational invariance of the static potential is 
observed. Using y/a = 440 MeV to set the scale, the lattice spacings for the 
coarse and fine lattices are 0.40(6) fm and 0.23(1) fm with x 2 /-/Vdf of 1.46 and 
0.91 respectively. We analyze 155 configurations on our coarse 6 3 x 12 lattice 
and 100 configurations on our fine 10 3 x 16 lattice. 

Five quark propagators are computed by the Stabilized Biconjugate Gra- 
dient algorithm EH for each configuration. The five quark masses selected are 
approximately 210, 180, 150, 120, 90 MeV, for both lattices. The second value, 
180 MeV, is taken as the strange quark mass. A point source is used at space- 
time location (x,y,z,t)=(l,l,l,2) on the 6 3 x 12 lattice and (1,1,1,3) on the 
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0.30 

Figure 3: Hadron masses in lattice units as a function of m q on our fine lattice. The lines 
are chiral fits as discussed in the text. For better viewing, the decuplet masses are shifted 
upward by 0.5 unit. 

10 3 x 16 lattice. The gauge-invariant smearing methodB is applied at the sink 
to increase the overlap of the interpolating operators with the ground states. 

Statistical errors are estimated in a third-order, single-elimination jack- 
knife, with bias corrections A third-order jackknife provides uncertainty 
estimates for the correlation functions, fits to the correlation functions, and 
quantities extrapolated to the chiral limit. 

3.2 Hadron Masses 

Fig. H shows the extracted hadron masses plotted as a function of the quark 
mass. We report values taken from covariance matrix fits to the time slice 
interval 4 through 8 on our coarse lattice and 6 through 9 on our fine lattice. 
These regimes provide the best signal-to-noise and good correlated xV-^df- 
n cr is determined by linearly extrapolating as a function of m q to zero. 
Similarly, the pseudoscalar kaon is extrapolated via m 2 K = cq + c\ m q . The 
form M = cq + c\ m q is used for all other extrapolations to the chiral limit. Fits 

3/2 

including an additional term C2m q are also considered and similar results are 
found with slightly larger error bars. 

Ratios of the chirally extrapolated masses are given in Table [| along with 
the ratios as observed in nature At = 7.0 the lattice spacing estimates 
follow the familiar pattern having the value based on the string tension lying 
between that of the p and nucleon based values. This is most likely an artifact 
of the quenched approximation. However, at (3 = 6.25 we find significant 
disagreement among the values and an unusual reordering of values. 

Focusing first on ratios of hadrons having the same angular momentum, 
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Table 1: Mass ratios after extrapolation to the chiral limit. a p and ajv are lattice spacings 
in fm set by the rho mass (770 MeV) and the nucleon mass (938 MeV). 



(3=6.25 


/3=7.0 


Expt. 


/3=6.25 


0=7.0 


Expt. 


Vector /Vector 






Pseudoscalar /Vector 




K*/p 1.20(2) 


1.20(2) 


1.16 


K/p 0.76(2) 


0.75(2) 


0.64 


<t>/p 1.40(3) 


1.39(3) 


1.32 








Octet/Octet 






Octet /Vector 






A/TV 1.16(2) 


1.16(1) 


1.19 


N/p 1.48(6) 


1.44(5) 


1.22 


E/7V 1.19(2) 


1.20(1) 


1.27 


k/p 1.72(6) 


1.66(5) 


1.45 


B/N 1.34(3) 


1.33(2) 


1.40 


£/p 1.76(6) 


1.72(5) 


1.55 








E/p 1.97(6) 


1.91(6) 


1.71 


Decuplet/Decuplet 






Decuplet /Vector 






S*/A 1.12(1) 


1.11(1) 


1.12 


A/p 1.79(7) 


1.77(5) 


1.60 


S*/A 1.24(3) 


1.22(2) 


1.24 


S*/p 2.00(6) 


1.97(6) 


1.80 


n/A 1.34(3) 


1.33(2) 


1.36 


H*/p 2.20(7) 


2.16(6) 


1.99 








fi/p 2.40(7) 


2.35(6) 


2.17 


String Tensions 






Decuplet / Octet 






a st 0.40(3) 


0.220(2) 




A/N 1.21(5) 


1.23(2) 


1.31 


a p 0.30(1) 


0.205(6) 




E*/7V 1.35(5) 


1.37(2) 


1.47 


ajv 0.36(1) 


0.242(6) 




B/N 1.49(4) 


1.50(3) 


1.63 








n/JV 1.63(4) 


1.64(3) 


1.78 



we see very little change in the values as the lattice spacing is decreased. These 
ratios are also remarkably similar to those observed in nature, despite the fact 
that these are quenched QCD calculations. In addition, these ratios support 
our selection for the strange quark mass. 

This close resemblance to nature is not shared by ratios of hadrons with 
different angular momentum. All four classes of ratios significantly disagree 
with those of nature. Once again we see the familiar quenched artifact of the 
Octet /Vector ratio being too large and the Decuplet/Octet ratio being too 
small. 

The standard failure of the Kj p mass ratio in the quenched approxima- 
tion is also seen here. This shortcoming has been widely realized through an 
examination of the J-parameterEj defined by 



dm p 

J = m„ ~—r 

ami 



m K » - m p 
~ m K ' — 2 T ■ ( 5 ) 

-i p /m x = 1.8 m K — TO^ 



Empirically this ratio is 0.48. However we find 0.42 on our coarse lattice and 
0.43 on our fine lattice. The physics associated with this discrepancy was 
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Figure 4: Plot of the p-meson mass as a function of the squared pion mass obtained from 
our finer lattice. a p has been used to set the scale. The dashed line illustrates the standard 
linear extrapolations of mfa and m p . The solid and dot-dash curves include the two- pion 
self-energy of the p mesonEd for dipole dispersion cut-off values of 1 and 2 GeV respectively. 
The increase in the slope at mp/m^ = 1.8 (m^ ~ 0.21) provided by the two-pion self energy 
is the right order of magnitude to restore agreement with the empirical value. 




first reported by Cohen and leinweber E3 where it was pointed out that the 
self-energy generated by two-pion intermediate states of the p-meson, which 
is excluded in the quenched approximation, acts to increase the J parameter. 
Fig. ^ provides a sketch of how including the two-pion self-energy of the p can 
increase the value of J to 0.46. 

Perhaps the most important information displayed in Table |l] is that the 
Octet/Vector mass ratios display less than satisfactory scaling for the larger 
lattice spacing. To further examine scaling and make contact with other stud- 
ies, we focus on the the N/ p mass ratio which is among the the most revealing 
of ratios. 

Fig. |B| shows a comparison of the N/ p mass ratio versus M p a at the chiral 
limit r-Fig- ^| shows the N/p mass ratio as a function of M p a at a fixed n/p mass 
ratioEa of 0.7. This method is free of complications from chiral extrapolations. 
Both cases clearly show the improvement provided by the D%34 action. Indeed, 
the D%34 action has reproduced the state-of-the-art quenched QCD ratios 
using unimproved actions at coarse lattice spacings of 0.24 fm. 



3.3 Dispersion and Rotational Symmetry 

In addition to mass ratios, hadron states at finite momentum projections pa = 
n(2n/L) are also calculated. Dispersion is examined by calculating the effective 
speed of light, defined by c 2 = (E 2 (p) — E 2 (0))/p 2 , which is to be compared 
with 1. _ 

A comparison with SW and D234 lattice actions B is made in Table 0. 
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Figure 5: The N/p mass ratio versus M p a at the chiral limit. Solid symbols denote t 
standard WHson action. Open symbols denote improved actions including SWLl (o), D234 
(A), Dx34t3 (□), and D*34 (o) (this work). 



1 



1.3 1 1 1 1 1 

0.0 0.5 1.0 1.0 2.0 

M a 

Figure 6: The N/p mass ratio versus M p a at a fixed tr/p mass ratio of 0.7 for various 
actions. The solid symbols denote the standard actions: Wilson (square and diamond), 
stagearcd (circle). The open symbols denote improved actions: nonperturbatively-improved 
Swla (x), fixed-point actionEil (□), SW (o), D234 (A), and D%34 (o) (this work). 



The dispersion for the 0(a)-improved SW action is very poor relative to the 
excellent dispersions of the next-nearest-neighbor improved Da;34 actions in 
general. The D%34 dispersion is excellent even at our coarse lattice spacing. 

Rotational symmetry is explored in Table [| At the coarser lattice spacing, 
some drift in the central values is seen for the pion arid nucleon. The drift in 
the pion is similar to that seen for the D234c jactionH. However, the drift in 
dispersion previously reported for the mesontj is not apparent in our results 
for the D^34 action- The D%34 action has much better rotational symmetry 



than the SW actional. The D%34 action provides satisfactory 
finer lattice spacing and is competitive with the D234 action! 



ispersion at our 
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Table 2: Comparison of SW, D234 and D%34 actions for the speed of light squared obtained 
from the dispersion of tt and p mesons at /m p ~ 0.7 for pa = (2ir/L). 



Hadron 


a (fm) 


SW 


D234 


D X 34 


TT 


0.40 


0.63(2) 


0.95(2) 


0.99(3) 


TT 


0.24 




0.99(4) 


1.04(4) 


P 


0.40 


0.48(4) 


0.93(3) 


0.93(6) 


P 


0.24 




1.00(6) 


0.99(6) 



Table 3: Evaluation of dispersion and rotational invariance via the effective speed of light. 
The results are for m q ~ 180 MeV. 



a (fm) 


n 


TT 


<t> 


N 


n 


0.40 


(1,0,0) 


0.98(2) 


0.91(4) 


1.00(7) 


0.99(12) 




(1,1,0) 


0.91(4) 


0.91(6) 


0.94(8) 


0.91(7) 




(1,1,1) 


0.86(9) 


0.92(10) 


0.92(7) 


0.90(10) 


0.24 


(1,0,0) 
(1,1,0) 

(1,1,1) 


1.04(3) 
1.05(4) 
0.98(6) 


1.02(4) 
1.02(4) 
0.98(6) 


1.10(6) 

1.06(5) 
1.06(6) 


1.06(7) 

1.11(5) 
1.06(5) 



4 Conclusion 

We have computed masses and dispersion relations of light hadrons in lat- 
tice QCD using tree-level 0(a 2 ) tadpole-improved gauge and fermion actions. 
These actions have the appeal of being simple to implement and inexpensive 
to simulate. A great deal of effort is being directed toward finding the ulti- 
mate improved action that will facilitate simulations on the coarsest of lattices. 
We note however, that many quantities of phenomenological interest such as 
hadron form factors involve momenta on the order of a GeV. As such, a highly 
improved action which is costly to simulate may not be the ideal action for 
hadron phenomenology, especially for exploratory purposes. 

The mass ratios obtained from the D%34 action at 0.24 fm on a modest 
10 3 x 16 lattice reproduce the state-of-the-art results using conventional unim- 
proved actions. Excellent dispersion and rotational invariance up to pa w 1 
are also found. These results demonstrate that the D%34 action can serve as 
a viable candidate for the study of hadron phenomenology, and in our view is 
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preferable to the highly-improved but more costly D234 action. We plan to 
use the D^34 action to study hadron properties beyond the spectrum, such as 
multipole form factors of hadrons in general. These results also bode well for 
future explorations beyond the quenched approximation. 

Support from the U.S. DOE and the Australian Research Council is grate- 
fully acknowledged. 
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